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Disruption of the high-affinity nitrate transporter NRT2.1 activates the priming defense 
against Pseudomonas syringae, resulting in enhanced resistance. In this study, it is 
demonstrated that the high-affinity ammonium transporter AMTl 1 is a negative regulator 
of Arabidopsis defense responses. The T-DNA knockout mutant amtl 7 displays enhanced 
resistance against Plectosphaerella cucumerina and reduced susceptibility to P syringae. 
The impairment of AMTll induces significant metabolic changes in the absence of 
challenge, suggesting that amtl.1 retains constitutive defense responses. Interestingly, 
amtll combats pathogens differently depending on the lifestyle of the pathogen. In 
addition, N starvation enhances the susceptibility of wild type plants and the mutant 
amtl 1 to P syringae whereas it has no effect on P cucumerina resistance. The metabolic 
changes of amtl 1 against P syringae are subtler and are restricted to the phenylpropanoid 
pathway, which correlates with its reduced susceptibility. By contrast, the amtl 1 mutant 
responds by activating higher levels of camalexin and callose against P cucumerina. In 
addition, amtll shows altered levels of aliphatic and indolic glucosinolates and other 
Trp-related compounds following infection by the necrotroph. These observations indicate 
that AMTl 1 may play additional roles that affect N uptake and plant immune responses. 

Keywords: AMT1.T, basal resistance, metabolomics, transceptor, NRT2.1 



INTRODUCTION 

Arabidopsis defense responses against P. syringae and P. 
cucumerina have been widely studied. These pathogens have very 
different life styles, in that P. syringae is a hemibiotrophic bacterial 
pathogen that penetrates through natural openings, whereas P. 
cucumerina is a necrotrophic fungus that can overcome cell barri- 
ers and penetrate into the cytoplasm (Berrocal-Lobo et al., 2002). 
Salicylic acid-dependent defenses have emerged as the most effec- 
tive resistance mechanism against the bacterium (Brooks et al., 
2005; Glazebrook, 2005). The effector secretion machinery of 
P. syringae is currently under study because of its effectiveness 
in suppressing Arabidopsis defenses. Although the direct link 
remains unknown, Camanes et al. (2012b) demonstrated that 
coronatine, a P. syringae effector, targets NRT2. 1 and subsequently 
causes a H2O2 burst. The nrt2 mutant remains insensitive to this 
manipulation and the increased JA-signaling that is responsive to 
coronatine (Melotto et al., 2006). P. cucumerina is more complex. 
Arabidopsis defenses against this pathogen are horizontal and 
involve several defense signaling pathways, such as SA, JA, lAA, 
and ABA- mediated callose deposition (Lorenzo and Solano, 2005; 
Garcia-Andrade et al., 2011; Sanchez- Vallet et al, 2012; Gamir 
et al, 2014). N fertilization influences plant-pathogen interac- 
tions; therefore, the following two factors affect plant resistance: 
the total amount of N and the source of N. Over-fertilization 
increases the severity of mildew in several grapevine cultivars and 
many crop plants (Keller et al., 2003; Marschner, 2012). The level 



of N in fertilization programs also affects plant resistance based 
on the pathogen lifestyle. For example, elevated levels of N pro- 
mote susceptibility against biotrophs, such as powdery mildew, 
whereas it results in reduced disease development of Botrytis 
cinerea (Walters and Bingham, 2007). However, the source of N 
appears to have a greater influence in plant resistance (Gupta 
et al., 2013). Fertilizing tobacco plants with NO3" accelerates 
the hypersensitive response following P. syringae pv. Phaseolicola 
infections and provides enhanced resistance. This resistance is 
related to NO3" reduction and NO production that triggers SA 
synthesis (Durner and Klessig, 1999). Plants fed NO3" display 
increased PRl expression and higher levels of SA. By contrast, 
tobacco plants fed with NH^ are more susceptible to the bac- 
terium. The bypass of nitrate reduction avoids NO production 
and triggers 4-aminobutyric accumulation, which is a nutrient 
for the pathogen (Solomon and Oliver, 2001), increasing the 
susceptibility of the plant. 

Previously, first in animal sciences and then in plant sciences, 
the term transceptor is applied to membrane proteins that per- 
form a dual transport and signaling function. Gojon et al. (2011) 
revised the transceptor role of the nitrate transporter family 
member NRTl. In addition to binding to low affinity nitrate, 
this protein can sense nitrate in the root environment. Additional 
studies demonstrated that both NRTl.l and the high-affinity 
transporter NRT2.1 are transceptors. The sensing network reg- 
ulated by NRTl.l has not yet been elucidated, although several 



www.frontiersin.org 



May 2014 | Volume 5 | Article 231 | 1 



Pastor et al. 



AMT1.1 displays enhanced pathogen resistance 



of its components have been recently described. One target of 
NRTl.l is the transcription of NRT2.1. The dual role of the 
induction of NRT2.1 during early time-points of nitrate sup- 
ply and repression during later time-points both directly and 
indirectly involve NRTl.l. Of the N metabolism components, 
mutants impaired in NRTl.l, such as chll are also affected in 
NIAl, NiR (encoding nitrate reductase and nitrite reductase, 
respectively) (Wang et al, 2009) and CIPK8, which has a reg- 
ulatory role on NRT2.1 (Hu et al, 2009). Glutamic acid (Glu) 
is important in N flow and also regulates root architecture by 
promoting root branching. This effect is antagonized by excess 
NO3", and this antagonistic action is also regulated by NRTl.l- 
dependent signaling (Walch-Liu and Forde, 2008). Root devel- 
opment and branching is also conditioned to aujdn signaling, 
and it appears that NRTl.l couples NO3" sensing and auxin 
signaling. 

There is recent evidence for other NO3" sensors coordinat- 
ing signaling events in Arabidopsis. NRT2.1 displays NRTl.l- 
independent roles and is also unrelated to NO3" uptake activity. 
This gene is a repressor of lateral root initiation under high 
sucrose and NO3" supply, but surprisingly, this gene seems to 
sense abiotic stress and concomitantly suppresses biotic stress 
responses. We have demonstrated that the mutant nrt2 is affected 
in bacterial effector manipulation (Camanes et al, 2012b). 
Mutation of NRT2.1 has significant consequences in the tran- 
scriptome of Arabidopsis upon P. syringae infection, suggesting 
a complex regulatory function for the NRT2. 1 gene and also the 
NRT2.1 protein (Camanes et al, 2012b). Therefore, we hypoth- 
esize that NRT2.1 should be designated a transceptor as previ- 
ously proposed by Gojon et al. (2011). In Arabidopsis, the link 
between the immune system and nitrate transport is not fuUy 
understood. However, recent advances indicate two different rel- 
evant events. First, low nitrate treatment induces a rapid burst 
of ethylene production and upregulates CTRl, EIN3, and EILl 
(Zheng et al, 2013). Second, we demonstrated that a low N treat- 
ment enhances P. syringae susceptibility (Camanes et al., 2012b), 
which may be related to the increase of ET-dependent signal- 
ing that is antagonistic to the SA-signaling that is recognized as 
an efficient defense against biotrophic pathogens (Pieterse et al., 
2009). 

Interactions with the plant immune system are not restricted 
to NRT2.1. Other NRT family members, such as NRT2.6 
show enhanced expression upon Erwinia amylovora infection 
(Dechorgnat et al, 2012). Despite high NO3" levels that trigger 
NRT2.6 expression, no nitrate-related phenotype has been asso- 
ciated with the nrt2.6 mutant. Furthermore, nrt2.6 is less tolerant 
to E. amylovora because of the reduced ability to produce H2O2 
upon infection. 

Some members of the AMTl family are high-affinity ammo- 
nium transporters. Of these, AMTl.l has the highest affinity for 
NH+ (Shelden et al., 2001). The removal of N increases AMTl.l 
expression, and this correlates with an increase in NH^ uptake. 
This function is affected by 30-40% in the insertional T-DNA 
mutant amtl.l. Interestingly, amtl.l plants are indistinguishable 
from wild type under optimal fertilization and growth condi- 
tions. Furthermore, amtl.l displays wild type levels of total N 
in the presence of N or after 4 days of N starvation (Kaiser 



et al, 2002). Although the ammonium supply may affect the 
plant responses to pathogens, to the best of our knowledge, 
no NH^ transporter is involved in the plant immune system. 
NRT2.1 and AMTl.l are the primary high-affinity transport sys- 
tems (HATS) involved in NO3" and NH^ uptake, respectively 
(Cerezo et al., 2001; Kaiser et al., 2002), and reciprocal regu- 
lation between AtNRT2.1 and AtAMTl.l expression has been 
observed (Camanes et al, 2012a). High-affinity ammonium and 
nitrate transporters act as sensors of low NH^ and NO3" in the 
root environment; therefore, thus the gene expression of AMTl.l 
and NRT2.1 is enhanced upon low fertilization. Interestingly, 
its role in HATS under normal N supply conditions ( > 1 mM 
NO3" or NH+) is negligible (Gansel et al., 2001). Using an 
experimental system with optimal fertilization, we focused on 
the role of the high-affinity ammonium transporter gene in the 
Arabidopsis immune responses against the hemibiotrophic bac- 
terium P. syringae and the necrotrophic fungus P. cucumerina. 
Metabolomic and genetic studies revealed an interplay between 
AMTl.l and plant resistance. 

MATERIALS AND METHODS 

PLANT MATERIAL AND GROWTH CONDITIONS 

Seeds of the Arabidopsis accession Col-0 were obtained from 
the SALK institute (Alonso et al, 2003), the EMS linl line (in 
Col-0 background) was generously supplied by Jocelyn Malamy 
(University of Chicago, EEUU, Little et al., 2005) and the line 
nrtS.l (saLk_043672) in the Col-0 background was obtained from 
SALK institute (Alonso et al, 2003). The background Col3-gll 
and the mutant amtl.l (in Col3-gll background) were provided 
by Alain Gojon (INRA, Montpelier, France). 

For the bioassays of bacterial resistance, 1 week after germina- 
tion, the seedlings were transferred to 33-ml soil pots. The plants 
were cultivated at 20°C day/18°C night with 8.5 h of light (105 [xE 
m^^ s^') per 24 h and 60% relative humidity. 

All plant genotypes were germinated in soU, and, 1 week 
after germination, seedlings were individually transferred to 
33 ml pots containing commercial potting soil (TKSl, Floragard 
GmbH; http://www.floragard.de). Plants were cultivated at 20°C 
day/ 1 8° C night temperatures with 8.5 h of light (105 |xE mT^ s^') 
per 24 h and 60% relative humidity. 

For fungal resistance assays. About 50 seeds were germinated 
in 33-ml soil pots. Plants were grown in the same conditions as 
described for bacterial experiments. 

PATHOGEN STRAINS AND RESISTANCE ASSAYS 

The bacterial strain Pseudomonas syringae pv. tomato DC3000 
was grown overnight at room temperature in King's B solid 
medium with appropriate antibiotics and diluted to the desired 
concentration with 10 mM MgS04 for plant inoculation. The 
bacteria were used to infect 5-week-old (or otherwise mentioned) 
Arabidopsis plants by dipping in a suspension of 2. 5 x 10^ colony- 
forming units (cfu)/ml using 0.02% Silwet L-77 (Tornero and 
Dangl, 2001), control plants were treated with Silwet accord- 
ingly. Three days after challenge inoculation, the disease level was 
determined by harvesting the plants a plating a range of dilu- 
tions in agar- KB medium containing rifampicin. After incubation 
at 28°C for 3 days, the number of rifampicin-resistant colony 
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forming units per gram of infected leaf tissue was determined, 
and bacterial proliferation over the 3-day time interval was 
calculated. 

The adapted strain of P. cucumerina BMM was kindly provided 
by Brigitte Mauch-Mani (Universite de Neuchatel, Switzerland). 
Two week old plants were sprayed with 1 x 10^ spores/mL and 
maintained at 100% relative humidity. Two days after inoculation 
leaves were sampled by freezing in liquid N2 and stored at — 80°C 
until analysis. Disease rate was determined by microscopical anal- 
ysis after tripan blue stainings of infected leaves (Flors et al., 
2008). Leaves were classified in a disease ranking as explained in 
the figure legend. 

GROWTH CONDITIONS FOR STARVATION EXPERIMENTS 

Arabidopsis thaliana was grown hydroponically as described in 
Lejay et al. (1999). The seeds were germinated directly on top of 
modified Eppendorf tubes filled with pre-wet sand. Plants were 
grown until the age of 5 weeks on a 1 mM NH4NO3 nutrient solu- 
tion (repressed plants), which prevented any growth difference 
between the two genotypes (Lejay et al., 1999). Before inoculation 
experiments, plants were transferred for 48 h to N-free solution 
(de-repressed plants). 

HPLC-ESI FULL SCAN MASS SPECTROMETRY (Q-TOF INSTRUMENT) 

Metabolome analysis was performed using an Acquity UPLC 
system (Waters, Mildford, MA, USA) interfaced to hybrid 
quadrupole time-of- flight (QTOF Premier). The LC separation 
was performed by HPLC SunFire C18 analytical column, 5 [im 
particle size, 2.1 x 100 mm (Waters). Analytes were eluted with 
a gradient of methanol and water containing 0.01% HCOOH. 
Chromatographic conditions and QTOF MS parameters were 
followed as described in Gamir et al. (2014). 

HPLC-ESI TANDEM MASS SPECTROMETRY (TRIPLE QUADRUPOLE 
INSTRUMENT). HORMONAL ANALYSIS 

An Acquity ultra-performance liquid chromatography system 
(UPLC) (Waters, Mildford, MA, USA) was interfaced to a triple 
quadrupole mass spectrometer (TQD, Waters, Manchester, UK). 
The chromatographic separation conditions were closely related 
to those described previously. Chromatograpic conditions and 
TQD parameters were followed as described in Flors et al. (2008) 
and Gamir et al. (2012). Masslynxv4.1 (Waters, Manchester, UK) 
software was used to process the quantitative data obtained from 
calibration standards and samples. 

FULL SCAN DATA ANALYSIS 

Raw data obtained from Masslynx software was transformed 
to.CDF using Databrigde provide by Masslynx package. The.CDF 
data was process with R for statistical computing using XCMS 
package for relative quantification (Smith et al., 2006). Principal 
Component Analysis (PGA) were used as explained in http:// 
www.numericaldynamics.com/ as a tool to define major changes 
in the metabolome of the plant under priming condition during 
fungal infections. 

For heatmap construction and clustering of metabolites it was 
used the software MarVis Filter and MarVis cluster (http://marvis. 
gobics.de/; Kaever et al., 2012). 



KINETICS OF 1^NH+ INFLUX 

The kinetics of ^^NH^ influx as a function of external ^^NH^ 
concentrations ([^^NH^]o) was measured with [^^NH^]o rang- 
ing from 0.02 to 0.8 mM. To kinetics studies control plants and 
de-repressed plants of three genotypes (Gol-0 and nrtS.l) were 
used. Influx of '^NH^ into the roots were assayed as described 
by Gazzarrini et al. (1999). Col-0 and nrt3.1 plants with a nor- 
mal fertilization or N depletion were sequentially transferred to 
0. 1 mM GaSQ4 for 1 min and to complete nutrient solution (pH 
6.0) containing '^NH;j" (98% atom excess ^^N) for 5 min, at the 
concentrations indicated in figures. At the end of the ^^N labeling, 
roots were washed for 1 min in 0.1 mM GaSQ4 and were sepa- 
rated from shoots. The roots were dried at 70° C for 48 h, weighed, 
crushed in a hammer-mUl and analyzed for total '^N content 
using an integrated system for continuous flow isotope ratio mass 
spectrometry (Euro-EA elemental analyser, EuroVector S.P.A.; 
and Isoprime mass spectrometer; GV Instruments). Root influx is 
expressed in jimol ^^NH^ (g root DW)^^ h^^. To kinetics stud- 
ies data-transformation method based on the Michaelis-Menten 
formalism was used. The experiment was repeated three times. 

RNA EXTRACTION AND REAL-TIME PCR ANALYSIS 

Gene expression by quantitative real-time RT-PGR was per- 
formed using RNA samples extracted from root tissue using 
the RNA kit (Omega Bio-Tek Inc, Doraville, GA, USA) accord- 
ing to the manufacturer instructions. To avoid contaminat- 
ing DNA, the samples were treated with DNAse I. A total of 
1 |xg of total RNA was annealed to oligo-dT and reverse tran- 
scribed using Omniscript Reverse Transcription kit (QIAGEN) 
to obtain cDNA. The sequences of the gene-specific oligonu- 
cleotides designed and used for real-time PCR are the following: 
AMTl.l forward: 5'acactgtggccagttaggcg3' and reverse: 5'ccgtg 
gggatgtctttgagaS', Tubuline (TUB) forward: 5'cgattccgttctcgatgt 
tgt3' and reverse: 5'aatgagtgacacacttggaatcctt3' and EFla forward: 
5'gtcgattctggaaagtcgacc3' and reverse: 5'aatgtcaatggtgataccacgc3'. 
Real-time PCR was conducted using the QuantiTect™ SYBR 
Green PCR Kit (QIAGEN) and the SmartCycler II instrument 
(Cepheid). The experiment was repeated three times. 

RESULTS 

ALTERATION IN NITRATE AND AMMONIUM TRANSPORTERS 
ENHANCES BASAL RESISTANCE 

It has been reported that deletions in NRT2.1 and NRT2.2 cause 
the reduced susceptibility of Arabidopsis to P. syringae (Camaiies 
et al, 2012b). NRT2.1 has been proposed as a transceptor that 
delivers and coordinates distal signaling that affects primary 
metabolism and defense signaling (Gojon et al., 2011; Camaiies 
et al., 2012b). The function of NRT2.1 in transporting nitrate 
is also influenced by the membrane protein NRT3.1, which is 
not a transporter itself, but its mutation significantly affects the 
HATS of nitrate (Okamoto et al, 2006; Orsel et al, 2006). To 
determine whether deficiencies in NRT3.1 and the ammonium 
transporter gene AMTl.l have consequences in the basal resis- 
tance, we characterized the responses of the mutant lines nrt3.1 
and amtl.l toward P. syringae and P. cucumerina. The enhanced 
resistance mutant linl (blocked in NRT2.1) was used as a con- 
trol. All three mutants showed enhanced resistance against both 
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P. syringae and P. cucumerina compared to their respective con- 
trol wild type plants (Figures 1A,B)- The mutant linl was more 
resistant to the bacterium, whereas nrt3.1 and amtl.l showed 
significantly reduced susceptibility, although lower than linl. 



Because NRT3. 1 influences nitrate transport, we tested whether 
it also affects ammonium transport. The mutant nrt3.1 showed 
wild type '^NHJ uptake kinetics both with N and after 48 h 
of N depletion (Figure 2A). In addition, we also determined 
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FIGURE 1 I Bacterial proliferation and disease rate in Col-0, linl, nrt3.1, 
and Col3-gl/1 and amtl.l plants infected with R syringae and R 
cucumerina respectivelly. (A) Five week old plants were 
challenge-inoculated by dipping in a bacterial suspension of R syringae at 
2 X 10^ c.f.u. mL"^ . The values presented are means (±SD) of the log of the 
proliferation values. Data represent the average of three independent 
experiments (n = 3). Different letters mean significant statistical differences 
(ANOVA, LSD test; p < 0.05). (B) Two week old plants were inoculated by 
spraying with 1 x 10-^ spores x mL^^ with PcBMM. Disease symptoms 



were recorded by trypan-blue staining at 5 days post inoculation. Disease 
rate was ranked according to the infected leaf surface: level I no infection, 
level II less than 25% of infected leaf surface, level III between 25 and 50% 
of infected leaf surface, level IV more than 50% of infected leaf surface. The 
figure shows a representative experiment that was repeated three times 
with the same results. Data presented are the means of the percentage of 
diseased leaves per plant. Asterisk indicate statistically significant differences 
compared with non-induced control plants (f-test; *p < 0.05; '^p < 0.05 with 
their respective controls; n~50 leaves). 
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FIGURE 2 I (A) Kinetics of ''^NH+ influx in Col-0 and nrt3. 1 roots in the low 
^^NH^ concentration range. Plants were grown hydroponically with N 
supplied as 1 mM NH4NO3 during 6 weeks. After that one group of plants 
were transferred to 1 mM NH4NO3, other group of plants were transferred to 
N-free nutrient solution (-N) during 48 h. ^^NHj" influx was measured at 
different concentrations of external ^^NH+. Each data is the mean of 30 
replicates ±SE. (B) Real-time PGR analysis of the expression of AMT1. 1 in 
Col-0 and nrt3. 1 plants fertilized normally along 5 weeks and exposed to N 
starvation 2 days (-N). The AMTl. 1 transcript levels were normalized to the 



expression of TUB measured in the same samples. The experiment was 
repeated using EFTa with similar results. Each bar represents average data 
with standard error bars from two technical replicates three independent 
experiments (n = 6). (C) Real-time PGR analysis of the expression of AMTl 7 
in mock and R syringae or R cucumerina wild type infected plants. The 
AMTl. 1 transcript levels were normalized to the expression of TUB measured 
in the same samples. The experiment was repeated using EFIa with similar 
results. Each bar represents average data with standard error bars from two 
technical replicates three independent experiments (n = 6). 



whether the AMTl.l gene is upregulated in the absence of N 
under the same experimental conditions. N depletion induces 
(de-represses) AMTl.l expression in wild type plants (Gazzarrini 
et al., 1999), whereas nrt3.1 shows wild type AMTl.l expression 
(Figure 2B). Both results suggest that NRT3.1 does not modulate 
ammonium kinetics uptake and that it does not alters AMTl.l 
de-repression. 

Finally, we tested whether infection with P. syringae or P. cuc- 
umerina affects AMTl.l expression in a manner that may act 
as targets for bacterial or fungal effectors. Following infection, 
AMTl.l expression was not significantly modified by any of the 
pathogens (Figure 2C). 

NITROGEN STARVATION ENHANCES P. syringae SUSCEPTIBILITY BUT 
DOES NOT AFFECT P. cucumerina RESISTANCE 

Previously, we determined that nitrate starvation for 48 h 
increased the susceptibility of Arabidopsis toward P. syringae and 



that nrt2.1 was also impaired (Camaries et al., 2012b). To deter- 
mine whether N depletion affects the basal resistance against both 
pathogens and to determine the contribution of the ammonium 
transporter AMTl.l to this interaction, we performed a hydro- 
ponic experiment by inducing N depletion for 48 h. According 
to our previous study (Camanes et al., 2012b), wild type starved 
plants were more susceptible to P. syringae (Figure 3), whereas 
they were not affected in their resistance to the fungus (Figure 3). 
Notably, amtl.l starved plants also showed increased susceptibil- 
ity against P. syringae, whereas their resistance to the fungus was 
not affected. 

amtl.l DISPLAYS ALTERED HORMONAL DEFENSE RESPONSES 
AGAINST P. syringae MIO P. cucumerina 

Because amtl.l showed reduced susceptibility and enhanced 
resistance toward the bacterium and the fungus, respectively, 
compared to the wild type plants, we analyzed the main 
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FIGURE 3 I R syringae proliferation and lesion diameter cause by R 
cucumerina in ColS-gll and amt7. 7 plants fertilized normally (Normal 
fertilization) along 5 weeks and exposed to N starvation 2 days before 
inoculation (N starvation). Hydroponically growing plants were 
challenge-inoculated either with a bacterial suspension of P syringae at 
2.5 X 10^ c.f.u./ml or with a drop of 1 x 10^ spores x mL"'' of PcBMM. 
Data represent the average of three independent experiments (n = 3). The 
values presented are means of infected plants (±SD). Asterisk indicates 
statistically significant differences (LSD test; p < 0.05). 



hormones involved in defense signaling. Surprisingly, the SA, JA, 
and JA-Ile levels remained lower in the amtl.l plants infected 
with P. syringae for the entire experiment compared to the wild 
type plants (Figure 4A). Indoleacetic acid (lAA), however, did 
not significantly change in the amtl.l plants. A reduction of JA 
and JA-Ile in amtl.l plants suggests that this mutation may coun- 
teract bacterial effectors by enhancing resistance. However, this 
hypothesis is unlikely because the SA levels remained lower in the 
mutant. Therefore, the hormonal analysis did not clearly explain 
the reduced susceptibility of amtl.l to the bacterium. By con- 
trast, the hormonal responses of the mutant against the fungus 
may contribute to resistance. JA-Ile and lAA following infection 
are increased compared to the wild type plants; however, SA at 
48 and 72 hpi is reduced in the mutant infected with P. cucume- 
rina (Figure 4B). We also tested other common defense responses 
that effective against P. cucumerina, such as the camalexin levels 
and callose accumulation (Ton and Mauch-Mani, 2004; Gamir 
et al., 2014). Both defensive responses were strongly enhanced in 
the mutant response to the pathogen. Camalexin and callose in 
amtl.l remained significantly higher at 48 hpi (Figure 5). 

METABOLOMIC PROFILING OF ainr7.nN RESPONSE TO PATHOGEN 
AHACK 

The hormonal profile of amtl.l in response to bacterial infection 
did not explain its enhanced resistance; therefore, we performed 



full metabolomic profiling to understand the metabolic changes 
that may contribute to the resistance of this mutant. For P. 
cucumerina, the classical defenses against necrotrophs, such as 
JA-Ile, callose, and camalexin, over-accumulated in amtl.l plants 
infected with the fungus. We also performed metabolic profiling 
following infection with the fungus to determine the extent to 
which amtl.l affects the responses to fungal infection. We used 
untargeted reverse-phase liquid chromatography-mass spectrom- 
etry (HPLC-QTOF-MS) for the profiling. CoB-gll and amtl.l 
plants were either mock or P. syringae or P. cucumerina infected. 
Metabolomic analysis and data reporting were performed as 
described by Pitzschke and Hirt (2010), Fernie et al. (2011), and 
Kaever et al. (2012). 

The acquired raw data were transformed into.cdf files using 
Databridge from the Masslynx 4.1 software (Masslynx 4.1, 
Waters). These data were subsequently subjected to analysis using 
the free software R for statistical purposes. The signals from the 
positive and negative electro-spray analysis (ESU-; ESI—) were 
processed separately, these ionization modes assure that all com- 
pounds that tend to form either cations or anions wiU be detected 
by the mass spectrometer. We performed a PCA of all signals 
obtained in ESU- and ESI—. The infection with P. syringae and 
P. cucumerina induced strong changes in the metabolome of the 
wild type plants and amtl.l plants that differed for the positively 
and negatively ionized compounds (Figure 6). The impact of the 
infection in ESU- appeared stronger in amtl.l upon P. cucume- 
rina infection compared with P. syringae; however the changes 
in ESI— are more subtle because there is an overlap in the vari- 
ability of amtl.l with or without infection independently of the 
pathogen. 

Notably, in the absence of infection, the basal metabolomes 
of CoB-gll and amtl.l do not overlap in any of the PCA 
with the exception of ESI— in the seedlings. This suggests that 
many metabolites are already altered in amtl.l in the absence 
of infection, and these changes may contribute to its enhanced 
resistance, discarding a general priming phenomenon to explain 
the phenotypes observed against both pathogens. Additionally, 
the wild type plants and amtl.l plants show differences in the 
basal responses between both experiments without inoculation, 
which may be due to a difference in the age of the plants 
(5 and 2 week-old plants for the P. syringae and P. cucume- 
rina experiments, respectively) but also to the different mock 
treatments (dip inoculated Silwet for P. syringae and sprayed 
10 mM MgS04 for P. cucumerina). A heatmap analysis was 
performed to identify metabolites involved in the resistance phe- 
notypes of amtl.l, and the clusters of metabolites that over- 
accumulated in amtl.l either in the absence or presence of the 
infection were selected (Figure SI). Clusters from the heatmaps 
that showed overaccumulated compounds in the mutant (red 
color) were selected for a detailed metabolic study and pathway 
analysis. 

Prior to identifying the signals in the selected clusters, we per- 
formed a detailed analysis of the amino acids in the metabolome. 
For such purpose, we constructed a library of amino acid stan- 
dards for the identity assignation using the exact mass provided 
by the Q-TOF analyser and the retention time parameter (Gamir 
et al., 2014). 
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FIGURE 4 I SA, JA, JA-lle, and lAA profiling upon P. syringae and R 
cucumerina infection. Plants were challenged as described in Figure 1. Both 
mock and pathogen infected plants were harvested at different time-points. 
Freeze dried material was processed for a targeted quantification analysis by 
TQD-MS. The concentration of the hormones was determined in all the samples 



by normalizing the chromatographic area for each compound with the dry 
weight of the corresponding sample. Leaf material from 15 individual plants for 
R syringae {A) resistance assays and 150 plants for /? cucumerina (B) resistance 
assays were pooled togetherfor each treatment x genotype combination. Data 
represent average three independent experiments ±SD; n = 3. 





CoB-gll 



Col3-gll 



amtl.1 



FIGURE 5 I Camalexin and callose levels upon /? syringae and R 
cucumerina infection. Plants were challenged as described in Figure 1 . Either 
mock or pathogen infected plants were harvested at 48 hpi. (A) Freeze dried 
material was processed for a targeted quantification analysis of camalexin by 
TQD-MS. The relative concentration was determined in all the samples by 
normalizing the chromatographic area for each compound with the dry weight of 
the corresponding sample. Leaf material from 1 50 individual plants were pooled 



together for each treatment x genotype combination. Data represent average 
three independent experiments. (B) Callose was visualized by aniline blue 
staining and epifluorescence microscopy (UV). Quantification was performed 
by determining the number of yellow pixels per million pixels corresponding to 
pathogen-induced callose on digital photographs of infected leaf areas. Asterisk 
indicates statistically significant differences (LSD test; p < 0.05). Data shown 
are means (±SD; n = 20) of the relative number of yellow pixels per photograph. 



In response to P. syringae, all the amino acids, together with 
azelaic acid and pipecolic, were not increased in the amtl.1 
mutant in response to the bacterium (data not shown). However, 
a different profile was observed following P. cucumerina infection. 



The majority of amino acids were increased in amtl.1 plants 
with or without infection (Figure 7). Pro, Thr, His, and Met were 
previously elevated in the absence of infection and this increase 
occurred until 48 hpi. Ala, Arg, Asn, Ile+Leu (not separated in 
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FIGURE 6 I Non-supervised Principal Component Analysis. (PCA) 
analysis representation of major sources of variability of ESI+ and ESI— 
signals obtained from a non-targeted analysis by HPLC-QTOFMS to 
monitor metabolomic changes during bacterial (A) and fungal invasion (B). 
(A) Five week old plants were dip inoculated with R syringae with 
2.5 X 10E^c.f.u/ml. 48hpi Leaf material from 15 individual plants were 
pooled together for each treatment x genotype combination. (B) Two 
week old plants were sprayed inoculated with lOE-^ spores/ml of R 



cucumerina and samples for analysis were collected 48 hpi. Leaf material 
from 150 individual plants were pooled together for each treatment x 
genotype combination. Data points represent two technical replicates 
from three independent experiments (biological replicates; n = 6) injected 
randomly into the HPLC-QTOFMS. The signals corresponding to different 
treatments were compared using the non-parametric Kruskal-Wallis test, 
and only data with a P-value lower than 0.01 between groups was used 
for subsequent processing. 



our chromatographic analysis), Cys, Gin, Lys, and Tyr showed 
a primed profile because they remained at the same levels in 
the wild-type plants before the inoculation; however, after infec- 
tion, amtl.l displayed elevated levels compared to Col3-gll. Only 
Asp, Glu, Phe, Trp, and Val in the mutant remained at the same 
levels as wild type infected plants (data not shown). Next, we 
performed a full comparative analysis of the metabolome. To 
identify the compounds, we used a library of standards with 
both the exact mass provided by the Q-TOF detector and the 
retention time. For the compounds with no available standards, 
we used the exact mass and fragmentation spectrum, when 
available, from the Massbank and Metlin databases. For such 
compounds, we assigned a tentative identification. After the sig- 
nals were either exactly or tentatively identified, the compounds 
were searched against the Aracyc (http://pathway.gramene.org/ 
gramene/aracyc.shtml) and Kegg (http://www.genome.jp/kegg/) 
databases to assign a putative metabolic pathway and a biolog- 
ical function. The analysis of the hormones or amino acids in 
amtl.l yielded a plausible explanation for its reduced suscep- 
tibility to P. syringae. However, the selection of the metabo- 
lite clusters from the heatmap that over- accumulated in amtl.l 
compared to the wild type plants following infection provides 
valuable information. Following the tentative identification of 
metabolites and their classification into metabolic pathways of 



Arabidopsis, we observed that fatty acid-CoA conjugates were 
increased in amtl.l either in the absence or presence of P. 
syringae (Figures). Nucleotides, such as ATP or UTP were also 
increased in amtl.l. Interestingly, three chemical species derived 
from TDP were also found for the amtl.l over- accumulated 
compounds. Dipeptides and tripeptides were another group of 
interest, which was surprising because the free amino acids in 
amtl.l are not different from the wild type plants. Although SA 
was not increased in the mutant, several benzoyl compounds 
over- accumulated in amtl.l, both with and without infection. 
The metabolic pathways with a higher number of hits in the clus- 
ter of over-accumulated compounds in amtl.l were the flavonoid 
and phenylpropanoid phytoalexins (Table 1). Although infec- 
tion induced a strong increase in these compounds in Col3-gll, 
amtl.l showed higher levels in the absence of infection. Finally, 
an aliphatic glucosinolate and vitamins B 1 and B6, all of which are 
involved in defense, were also increased in the mutant (Figure 8). 
For the interaction with P. cucumerina, the identified compounds 
matched four secondary metabolism pathways of Arabidopsis, 
indolic glucosinolates, aliphatic glucosinolates, the Trp pathway 
and phenylpropanoids (with a total of 19 compounds) (Figure 9). 
Notably, many secondary metabolites in these groups accumu- 
lated in amtl.l without infection, and all were then either 
increased or maintained in response to P. cucumerina in amtl.l 
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FIGURE 7 I Aminoacid profiling upon R cucumerina infection. Two week 
old Col3-gl/1 and amtl. 7 plants either mock or R cucumerina inoculated were 
processed for relative quantification analysis by HPLC-QTOFMS data. The 
concentration of the metabolites was determined in all the samples by 
normalizing the chromatographic area for each compound with the dry 



weight of the corresponding sample. White bars are mock inoculated and 
filled bars are P cucumerina infected plants. Leaf material from 150 individual 
plants were pooled together for each treatment x genotype combination. 
Boxplots represent average three independent experiments with two 
technical replicates (n = 6). 



plants compared to wild type plants. In addition to these path- 
ways, although there were fewer hits, tentative identification of 
the metabolites showed that amtl.l responded to infection by 
increasing three mevalonic acid derivatives, four nucleotide inter- 
mediaries, three shikimic acid metabolites, two polyamines, and 
several hormones, such as GA3, GA-A58, and ABA (Figure S2). 

DISCUSSION 

We previously reported (Camafies et al., 2012b) that under opti- 
mal nutritional conditions, the total contents of N in nrt2.1 
and amtl.l plants are identical, which indicates that the N con- 
tent is not responsible for altering the basal resistance in the 
mutants. Based on our previous study, we performed all resis- 
tance experiments, except otherwise noted, under the optimal 
fertilization conditions. First, we demonstrated that mutations in 
several genes related to N transport resulted in an increased resis- 
tance against two pathogens with different lifestyles. The HATS 
of nitrate uses a membrane complex in which the NRT3.1 pro- 
tein modulates the function of NRT2. 1 , and mutation of NRT3. 1 



results in the most significant losses of HATS function compared 
to NRT2.1 mutations (Orsel et al., 2006). Therefore, we tested two 
different hypotheses. First, we tested whether this gene may also 
regulate ammonium transport, and second, we tested whether 
NRT3.1 is involved in the plant defense response. The uptake 
experiments with '^NH^ showed that the absorption of ammo- 
nium in nrtS.l does not differ from Col-0 plants, which suggests 
that it is not involved in ammonium HATS. We also demonstrated 
that AMTl.l gene expression with or without N was similar for 
Col-0 and nrtS.l plants. These observations confirm that this 
gene is not involved in AMTl.l modulation. Because nrt3.1 plants 
show resistance against P. cucumerina and R syringae, this sug- 
gests that it may have overlapping mechanisms with NRT2.1 to 
stimulate defense; however, this hypothesis must be confirmed. 

We ascertained the influence of AMTl.l on plant resistance 
because our experiments suggest that AMTl.l acts as a tran- 
sceptor to coordinate plant tolerance against N depletion and 
defense against pathogens as suggested previously for NRTl 
and NRT2.1 (Gojon et al, 2011; Camafies et al, 2012b). The 
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FIGURE 8 I Profiling of the main hits in the overaccumulated 
compounds upon P. syringae (Pst) infection in amtl.1 compared with 
Col3-gl1. Five week old Col3-gl1 and amtll were infected as described 
in Figure 1. After 48hpi plants were processed for relative quantification 
analysis by HPLC-QTOFMS data. The concentration of the metabolites 
was determined in all the samples by normalizing the chromatographic 
area for each compound with the dry weight of the corresponding 
sample. The compounds were tentatively identified using the exact mass 



criteria using the METLIN and Massbank databases. The compounds 
were grouped by metabolic pathways according to KEGG and AraCyc 
databases. White bars are mock inoculated and filled bars are P syringae 
infected plants. Leaf material from 15 individual plants were pooled 
together for each treatment x genotype combination. Boxplots represent 
average three independent experiments with two technical replicates. 
Only data showing a p-value below 0.05 after a Kuskal Wallys test were 
used for pathway assignation (n = 6). 
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FIGURE 9 I Profile of the main hits in the overaccumulated 
compounds upon R cucumerina infection in amtl.l compared with 

Col3-gl1. Two week old Col3-gl1 and amtll were infected as described 
in Figure 1. After 48hpi plants were processed for relative quantification 
analysis by HPLC-QTOFMS data. The concentration of the metabolites 
was determined in all the samples by normalizing the chromatographic 
area for each compound with the dry weight of the corresponding 
sample. The compounds were tentatively identified using the exact mass 
criteria using the IVIETLIN and IVIassbanl< databases. Those compounds 



tentatively identified have the exact mass indicated. Te compounds fully 
identified have no reference to the mass. The compounds were grouped 
by metabolic pathways according to KEGG and AraCyc databases. White 
bars are mock inoculated and filled bars are P cucmerina infected plants. 
Leaf material from 150 individual plants were pooled together for each 
treatment x genotype combination. Boxplots represent average three 
independent experiments with two technical replicates. + indicates 
outliers. Only data showing a p-value below 0.01 after a Kuskal wallys 
test were used for pathway assignation (n = 6). 
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mutant amtl.l showed reduced susceptibility against P. syringae, 
and surprisingly, amtl.l showed strong resistance toward P. cuc- 
umerina infection compared to wild type plants. Previously, 
we demonstrated that environmental factors contribute to the 
enhanced expression of NRT2.1. For example, N depletion affects 
the resistance of Arabidopsis to P. syringae. Although the deple- 
tion of N induces AMTl.l gene expression (Figure!), it also 
increased the susceptibility in amtl.l plants, suggesting that 
it is the absence and not the overexpression of AMTl.l that 
alters resistance mechanisms. Notably, neither infection with P. 
syringae or P. cucumerina significantly altered the AMTl.l tran- 
scripts. We tested whether this also applies to resistance against 
necrotrophs. The enhanced resistance of amtl.l against P. cuc- 
umerina was not altered by N depletion; therefore, the absence of 
these genes impacts Arabidopsis defense against this pathogen but 
not their overexpression, suggesting posttranscriptional or post- 
translational regulation. This is not so surprising because NRTl.l 
and NRT2.1 have a complex protein regulation mechanism medi- 
ated by proteins, such as CIPK23 and NAR2.1 (Okamoto et al., 
2006; Ho et al., 2009). 

To identify the molecular mechanisms for the amtl.l 
phenotypes, we performed several targeted and non-targeted 
metabolomic studies. For the targeted studies, hormones reg- 
ulate the plant immune system (Pieterse et al, 2009). Active 
SA-dependent responses against P. syringae promote resistance 
in Arabidopsis. Additionally, because of the existing negative 
crosstalk between the SA and JA signaling pathways, resistance 
against biotrophs is associated with the concomitant downreg- 
ulation of JA and JA-Ile levels (Pieterse et al, 2009). This has 
been observed in amtl.l plants upon P. syringae infection. In fact, 
SA, JA, and JA-Ile remain lower compared to wild type plants 
during a time-course experiment. Therefore, the reduced suscep- 
tibility of amtl.l to the bacterium cannot be caused by changes 
in the hormonal balance. Except for lAA, the hormones tested 
remained at lower levels in the absence of infection. This suggests 
that disruption of AMTl.l affects the basal levels of SA, JA, and 
JA-Ile, producing several constitutive changes. Interestingly, these 
basal changes are not visible in 2-week old plants; therefore, some 
alterations affected by the mutation are age dependent. For non- 
host resistance, such as for P. cucumerina, the defense responses 
are horizontal, multigenic and not controlled by single hormonal 
pathway, because other hormones in addition to SA, JA, ABA, 
and ET influence Arabidopsis responses against this necrotroph 
(Sanchez-Vallet et al, 2010). These compounds include the P- 
subunit of the heterotrimeric G-protein (Delgado-Cerezo et al., 
2012) and glucosinolates and other Trp derivatives (Sanchez- 
Vallet et al., 2012; Gamir et al., 2014). In our first approach, we 
showed that amtl.l influences the SA, JA-Ile, and lAA profile 
during infection. Both JA-Ile and lAA over-accumulate, whereas 
SA is at lower levels compared to Col3-gll. This may partially 
explain the observed phenotype because P. cucumerina has a 
necrotrophic lifestyle and is sensitive to JA-dependent signaling 
(Thomma et al., 2000). Recent research has suggested complex 
interplay between hormones in addition to the negative crosstalk 
between SA and JA (Pieterse et al, 2009), lAA also has antago- 
nistic effects with SA, and both may explain the lower amounts 
of SA observed in amtl.l, although the final link between the 



ammonium transporter and hormonal signaling remains to be 
elucidated. Arabidopsis can also resist P. cucumerina through two 
other major mechanisms, caUose deposition and camalexin (Ton 
and Mauch-Mani, 2004; Garcia-Andrade et al., 2011; Sanchez- 
Vallet et al., 2012). Both mechanisms are directly or indirectly 
influenced by AMTl.l. Early callose deposition is enhanced in 
the mutant, and the camalexin levels remain higher at 48hpi 
compared to Col3-gll. 

In oui Arabidopsis- P. syringae system, we observed that amtl.l 
does not influence the basal amino acid accumulation in response 
to the bacterium. By contrast, there is a significant change in 
the amino acid profile for P. cucumerina. The mutant amtl.l 
over-accumulates several amino acids in the absence of infection, 
but these changes are even more robust after infection. Up to 
11 amino acids over-accumulate in amtl.l infected plants com- 
pared to Col3-gll. Surprisingly, Trp remains at the same level as 
the wild type plants with or without P. cucumerina. This may be 
explained by the over-accumulation of other Trp derivatives, such 
as indolic glucosinolates, lAA, camalexin, and other indole conju- 
gates (Figure 9). These compounds are metabolic sinks that may 
trap the putative overproduction of Trp. 

According to the PGA, the distance between Col3-gll and 
amtl.l upon infection by both pathogens is rather marked; 
however, the diversity of the metabolic responses against P. cuc- 
umerina is much wider compared with the bacterium. In fact, 
of all of the compounds analyzed by targeted and untargeted 
chromatography that were found in clusters of over-accumulated 
metabolites in amtl.l upon P. syringae infection, only five path- 
ways were identified, and of 45 signals, 24 were tentatively 
assigned to the phenylpropanoids (Table 1) and 5 to benzoic acid 
derivatives. This suggests that shikimic acid derivatives may com- 
pete with SA accumulation, which is a likely explanation for the 
reduced susceptibility of amtl.l and the reduced levels of free 
SA. Phenylpropanoids have been associated with plant resistance 
against a wide range of pathogens (Shadle et al., 2003). Two 
other compounds over- accumulated in amtl.l plants are vita- 
mins Bl and B6. Both have been linked to disease resistance 
(Denslow et al., 2005; Ahn et al., 2007). Exogenous treatment with 
vitamin Bl can induce defense priming in Arabidopsis against 
P. syringae mediated by H2O2 accumulation, callose deposition, 
and NPRl. Vitamin B6 is involved in superoxide quenching and 
stress responses, and during P. syringae infection in Nicotiana, 
vitamin B6 acts as an antioxidant and modulator of active oxy- 
gen species (Denslow et al., 2005). The relevance of the other 
metabolites, such as small peptides, nucleotides and GoA con- 
jugates in the defense against pathogens must be clarified. Of 
the signals over- accumulated in amtl.l infected with P. cucume- 
rina, 51 were tentatively or fuUy identified. These compounds 
were mainly involved in the indolic and aliphatic glucosino- 
lates pathways, Trp pathway, flavonoids, nucleotides, and fatty 
acids. Notably, glucosinolates have been directly linked to resis- 
tance against P. cucumerina Sanchez-Vallet et al., 2012; Gamir 
et al., 2014, and also callose deposition in response to PAMPs 
(Glay et al., 2009). Although no direct evidence for AMTl.l 
and glucosinolates has been previously reported, this may be 
worth investigating. The NRT/PTR transporters have been asso- 
ciated with glucosinolate translocation in Arabidopsis, mainly 
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aliphatic (Nour-Eldin et al, 2012). Additionally, Met is a precur- 
sor of aliphatic glucosinolates (Kraker and Gershenzon, 2011). 
This amino acid highly accumulates in amtl.l plants infected 
with P. cucumerina. Furthermore, three precursors of aliphatic 
glucosinolates are also elevated in amtl.l (Figure 9). These find- 
ings suggest a repressive influence of AMTl.l on plant defenses 
against necrotrophic infection that result in changes in amino 
acid metabolism and/or peptide transport that causes increased 
defenses, most likely mediated by aliphatic and indolic glucosi- 
nolates. Other relevant metabolic contributions from different 
pathways may also be involved, such as chlorogenic acid, ABA, 
spermidine, several fatty acids, and other unclassified compounds 
that are also over- accumulated in amtl.l plants upon infection. 

We previously demonstrated that the member of the high- 
affinity transport system of nitrate, the NRT2.1 gene, in addi- 
tion to transporting nitrate, is involved in the coordination of 
plant defense responses (Camanes et al., 2012b). This gene is 
de-regulated upon low nitrate concentrations in the roots, sub- 
sequently, the HATS is activated, providing tolerance mechanism 
against abiotic stresses (Gansel et al., 2001). Despite these func- 
tions, NRT2.1 represses plant responses against bacterial and 
fungal pathogens. The metabolic interplay between nitrate uptake 
and resistance against pathogens is being elucidated. Because this 
resistance mechanism is multicomponent and non-pathogen- 
specific, the disruption oi NRT2.1 results in enhanced resistance 
against R syringae and P. cucumerina through different defense 
mechanisms (Camanes et al., 2012b; Gamir et al., 2014). 

Here, we demonstrate that AMTi.i acts as a negative regulator 
of Arabidopsis defense. Its disruption has a moderate impact on 
metabolomic changes upon P. syringae infection, which is corre- 
lated with a reduced susceptibility to the bacterium. By contrast, 
the metabolomic changes of amtl.l upon P. cucumerina chal- 
lenge are significant, targeting Trp-derivatives and indolic glu- 
cosinolates. These changes may be responsible for the enhanced 
resistance of amtl.l against the fungus. This study and previ- 
ous findings from our laboratory indicate that N transporters 
are transceptors that impact the transcriptome (Camanes et al., 
2012b) and metabolome oi Arabidopsis (Gamir et al., 2012, 2014) 
to coordinate abiotic stress tolerance and the negative regula- 
tion of biotic stress responses. The manner in which NRT2.1 and 
AMTl.l regulate stress signaling remains to be identified; how- 
ever, it is clear that the plant has mechanisms to distinguish and 
regulate different sources of stress. Transceptors appear to play 
a key role in integrating the environmental signals and plant 
defense responses. 
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Figure SI | Heat map analysis performed withi IVIarvis (Filter and Cluster 
pacl<ages). Five and two week old plants either mock {Col3-gl/1 and 
amtl. 1) or either R syringae Col3-gl/1 P. syringae and amtt 1 R syringae or 
R cucumerina inoculated {Col3-gi/1 Rc and amtl 1 Rc) plants were 
processed for relative quantification analysis by HPLC-QTOFMS data. The 
concentration of the metabolites was determined in all the samples by 
normalizing the chromatographic area for each compound with the dry 
weight of the corresponding sample. Heatmaps are generated by using 
Mar-Vis Filter and Cluster following a Kruskal-Wallys test (p < 0.01 ). 
Clusters overepresented amtl. 1 mock or infected (intense red colors) 
compared with wild type plants were selected for subsequent data 
analysis. 

Figure S2 | Profile of additional hits in the overaccumulated compounds 
upon R cucumerina infection in amtl. 1 compared with Col3-gl1. For 

experimental details see Figures 1, 9. The concentration of the 
metabolites was determined in all the samples by normalizing the 
chromatographic area for each compound with the dry weight of the 
corresponding sample. The compounds were tentatively identified using 
the exact mass criteria using the METLIN and Massbank databases. 
Those compounds tentatively identified have the exact mass indicated. Te 
compounds fully identified have no reference to the mass. The 
compounds were grouped by metabolic pathways according to KEGG and 
AraCyc databases. White bars are mock inoculated and filled bars are R 
cucmerina infected plants. Leaf material from 150 individual plants were 
pooled together for each treatment x genotype combination. Boxplots 
represent average three independent experiments with two technical 
replicates. Only data showing a p-value below 0.01 after a Kuskal wallys 
test were used for pathway assignation (n = 6). 
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